TOPICAL REVIEW 



Wigner crystal physics in quantum wires 

Julia S Meyer 

Department of Physics, The Ohio State University, Columbus, Ohio 43210, USA 
K A Matveev 

Materials Science Division, Argonne National Laboratory, Argonne, Illinois 60439, 
USA 

E-mail: jmeyer@mps . ohio-state . edu 

Abstract. The physics of interacting quantum wires has attracted a lot of attention 
recently. When the density of electrons in the wire is very low, the strong repulsion 
between electrons leads to the formation of a Wigner crystal. We review the rich 
spin and orbital properties of the Wigner crystal, both in the one-dimensional and 
quasi-one-dimensional regime. In the one-dimensional Wigner crystal the electron 
spins form an antiferromagnetic Heisenberg chain with exponentially small exchange 
coupling. In the presence of leads the resulting inhomogeneity of the electron density 
causes a violation of spin-charge separation. As a consequence the spin degrees of 
freedom affect the conductance of the wire. Upon increasing the electron density, the 
Wigner crystal starts deviating from the strictly one-dimensional geometry, forming 
a zigzag structure instead. Spin interactions in this regime are dominated by ring 
exchanges, and the phase diagram of the resulting zigzag spin chain has a number 
of unpolarized phases as well as regions of complete and partial spin polarization. 
Finally we address the orbital properties in the vicinity of the transition from a one- 
dimensional to a quasi-one-dimensional state. Due to the locking between chains in 
the zigzag Wigner crystal, only one gapless mode exists. Manifestations of Wigner 
crystal physics at weak interactions are explored by studying the fate of the additional 
gapped low-energy mode as a function of interaction strength. 
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1. Introduction 

First experiments [U [2] on electronic transport in one-dimensional conductors revealed 
the remarkable quantization of conductance in multiples of the universal quantum 
2e 2 /h, where e is the elementary charge and h is Planck's constant. These experiments 
were performed by confining two-dimensional electrons in GaAs heterostructures to one 
dimension by applying a negative voltage to two gates, thereby forcing the electrons to 
flow from one side of the sample to the other via a very narrow channel. Such devices, 
typically referred to as quantum point contacts, are the simplest physical realization of 
a one-dimensional electron system. Although the length of the one-dimensional region 
in quantum point contacts is relatively short, the quantization of conductance indicates 
that transport in such devices is essentially one-dimensional. Longer quantum wires 
have been created later using either a different gate geometry [3], or by confining two- 
dimensional electrons by other means, such as in cleaved-edge-overgrowth devices [I]. 
Finally, a fundamentally different way of confining electrons to one dimension has been 
recently realized in carbon nanotubes [HIE]. The interest in the study of one- dimensional 
conductors is stimulated by the relatively low disorder in these systems and by the ability 
to control their parameters. For instance, the effective strength of the electron-electron 
interactions is determined by the electron density, which can be tuned by changing 
the gate voltage. Thus quantum wire devices represent one of the simplest interacting 
electron systems in which a detailed study of transport properties can be performed. 

Interactions between one-dimensional electrons are of fundamental importance. 
Unlike in higher- dimensional systems, in one dimension the low-energy properties of 
interacting electron systems are not described by Fermi-liquid theory. Instead, the 
so-called Tomonaga-Luttinger liquid emerges as the proper description of the system 
in which, instead of fermionic quasiparticles, the elementary excitations are bosons 
[Tj. Interestingly, the quantization of conductance in quantum point contacts is well 
understood in the framework of noninteracting electrons [E] despite the relatively strong 
interactions in these devices. This paradox was resolved theoretically [91 [10l [11] by 
considering a Luttinger liquid with position-dependent parameters chosen in a way that 
models strongly interacting electrons in the quantum wire connected to leads in which 
interactions can be neglected. It was found that the dc conductance of such a system is 
completely controlled by the leads, and is therefore insensitive to the interactions. 

The latter conclusion is in apparent disagreement with experiments observing the 
so-called 0. 7 structure in the conductance of quantum point contacts [121 E2 ELU EEl EEl 
rm [T8[ [191 [20] . This feature appears as a quasi-plateau of conductance at about 0.7 x 
2e 2 /h at very low electron density in the wire, and usually grows with temperature. A 
number of possible explanations have been proposed, most of which attribute the feature 
to the fact that at low densities the effective interaction strength is strongly enhanced. 
One of the most common explanations attributes the 0.7 structure to spontaneous 
polarization of electron spins in the wire [121 1131 [131 [IZl [ISl 1201 12H 12^1 1231 1231 12SI - 
Although such polarization is forbidden in one dimension [26], the electrons in quantum 
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Figure 1. (a) A onc-dimcnsional Wigner crystal formed in a quantum wire at low 
electron density, (b) The zigzag Wigner crystal forms in a certain regime of densities 
when the electrons are confined to the wire by a shallow potential. 



wires are, of course, three-dimensional, albeit confined to a channel of small width. This 
deviation from true one-dimensionality may, in principle, give rise to a spin-polarized 
ground state of the interacting electron system. 

The electrons in quantum wires interact via repulsive Coulomb forces. As a result of 
the long-range nature of the repulsion, at low density the kinetic energy of the electrons 
is small compared to the interactions. To minimize their repulsion, electrons form 
a periodic structure called the Wigner crystal [27]. In one dimension the long-range 
order in the Wigner crystal (Figure [Qa)) is smeared by quantum fluctuations [28], and 
therefore the crystalline state can be viewed as the strongly-interacting regime of the 
Luttinger liquid. However, the presence of strong short-range order provides a clear 
physical picture of the strongly interacting one- dimensional system and enables one to 
develop a theoretical description of quantum wires in the low-density regime. 

In the Wigner crystal regime the electrons are strongly confined to the vicinity of the 
lattice sites. As a result the exchange of electron spins is strongly suppressed, and only 
the nearest neighbor spins are coupled to each other. One can then think of the electron 
spins forming a Heisenberg spin chain with a coupling constant J much smaller than 
the Fermi energy Ep. The presence of two very different energy scales Ep and J for the 
charge and spin excitations distinguishes the strongly interacting Wigner crystal regime 
from a generic one-dimensional electron system with moderately strong interactions. In 
particular, the Luttinger liquid theory is applicable to the Wigner crystal only at the 
lowest energies, e < J. On the other hand, if any of the important energy scales of 
the problem exceed J, the spin excitations can no longer be treated as bosons, and the 
conventional Tomonaga-Luttinger picture fails. One of the most interesting examples of 
such behavior occurs when the temperature T is in the range J < T <C Ep. In this case 
the charge excitations retain their bosonic properties consistent with Luttinger liquid 
theory, whereas the correlations of electron spins are completely destroyed by thermal 
fluctuations. Such one-dimensional systems are not limited to the Wigner crystal regime 
and are generically referred to as spin-incoherent Luttinger liquids. We argue in Sec. [2] 
that the coupling of spin and charge excitations in this regime leads to a reduction of the 
conductance of the quantum wire from 2e 2 / h to e 2 / h. A number of additional interesting 
properties of spin- incoherent Luttinger liquids are discussed in a recent review [29|. 

The electrons in a quantum wire are confined to one dimension by an external 
potential. In the common case of the potential created by negatively charged gates 
placed on top of a two-dimensional electron system, the confining potential can be 
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rather shallow. In this case the strong repulsion between electrons can force them 
to move away from the center of the wire, transforming the one-dimensional Wigner 
crystal to a quasi-one- dimensional zigzag structure, Figure Hfb). In the case of classical 
electrons such a transition has been studied in [301 ED E2]; we review this theory in 
Sec. [3j The zigzag Wigner crystal has rich spin properties due to the fact that each 
electron can now be surrounded by four neighbors with significant spin coupling. Ring 
exchange processes play an important role and may under certain circumstances give 
rise to a spontaneous polarization of electron spins. The spin properties of the zigzag 
Wigner crystals are discussed in Sec. HI 

The transformation of a one-dimensional Wigner crystal to the zigzag shape is a 
special case of a transition from a one-dimensional to a quasi-one dimensional state of 
electrons in a quantum wire. Another such transition occurs in the case of noninteracting 
electrons when the density is increased until population of the second subband of 
electronic states in the confining potential begins. These two transitions seem to have 
rather different properties. Indeed, in the case of noninteracting electrons the population 
of the second subband entails the emergence of a second acoustic excitation branch in 
the system. On the other hand, even though the zigzag crystal has two rows, their 
relative motion is locked, and one expects to find only one acoustic branch in this case. 
It is therefore interesting to explore how the number of acoustic excitation branches 
changes as the interaction strength is tuned. In the regime of strong interactions this 
requires developing the quantum theory of the transition from a one-dimensional to 
a zigzag Wigner crystal. We discuss such a theory in Sec. [5j where it is shown that 
quantum fluctuations do not lead to the emergence of a second acoustic branch in the 
zigzag crystal. This feature of the Wigner crystal survives even at weak interactions, 
with the second acoustic branch appearing only when the interactions are completely 
turned off. 

2. One-dimensional crystal 

2.1. Quantum wire at low electron density 

Electrons in a quantum wire repel each other with Coulomb forces. To characterize 
the strength of interactions, let us compare the typical kinetic energy of an electron, 
which is of the order of the Fermi energy Ep ~ h 2 n 2 /m, with the typical interaction 
energy e 2 n/e. (Here n is the electron density, m is the effective mass, and e is the 
dielectric constant of the medium.) Clearly, the Coulomb repulsion dominates over the 
kinetic energy in the low-density regime nas <C 1, where a# = h 2 e/me 2 is the Bohr's 
radius in the material. Then the ground state of the system is achieved by placing 
electrons at well-defined points in the wire, separated from each other by the distance 
n -1 , Figure [H(a), thus creating a Wigner crystal. Because the kinetic energy of electrons 
is small, the amplitude 5x of the zero-point fluctuations of electrons near the sites of 
the Wigner lattice is much smaller than the period of the crystal, n5x ~ (na B ) 1//4 <C 1. 
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In experiment the quantum wire is usually surrounded by metal gates. As a result, 
the Coulomb interactions between electrons are screened at large distances by image 
charges in the gates. For example, if the gate is modeled by a metal plane at distance 
d from the wire, the interaction potential becomes 



At large distances this potential falls off as V(x) ~ 2e 2 d 2 /e\x\ 3 , much more rapidly than 
the original Coulomb repulsion. As a result, in the limit n —>■ the crystalline ordering 
of electrons will be destroyed by quantum fluctuations. Comparison of the Fermi energy 
with the screened Coulomb repulsion ([1]) shows that the Wigner crystal exists in the 
range of densities asd~ 2 <n< a^ 1 , provided that the distance to the gate d^> a B . In 
typical experiments with GaAs quantum wire devices as = lOnm and d > lOOnm; thus 
the Wigner crystal state should persist until unrealistically low densities ~ 10 _3 nm _1 . 

Similar to phonons in conventional crystals, the Wigner crystal supports acoustic 
plasmon excitations — propagating waves of electron density. The speed of plasmons is 
given by[j] 



The Hamiltonian describing these low-energy excitations is easily obtained by treating 
the Wigner crystal as a continuous medium. Adding the kinetic energy and the potential 
energy of elastic deformation, one obtains 



where u(x) is the displacement of the medium at point x from its equilibrium position, 
and p(x) is the momentum density. In one dimension the acoustic excitations destroy 
the long range order in the crystal even at zero temperature, {[u(x) — u(0)] 2 ) ~ 
(h/irmns) lnnxj§| 

In the model of spinless electrons, Hamiltonian ([3]) accounts for all possible low- 
energy excitations of the system. However, in the presence of spins, there are additional 
excitations not included in In the Wigner crystal regime the electrons are localized 
near their lattice sites, Figure [H(a), and to a first approximation the spins at different 
sites are not coupled. The exchange coupling of two spins at neighboring sites occurs 
via the process of two electrons switching their places on the Wigner lattice. When the 
electrons approach each other, the strong Coulomb repulsion creates a high potential 
barrier. As a result, the exchange processes are very weak, and only the coupling of the 

| The result ((2]) was derived in [33] for densities in the range dT 1 n <C a^. Extending their 
calculation to the density range asjd 2 <§; n <§; dT 1 , one finds s = [24e 2 n 3 c? 2 C(3)/em] 1 / 2 . 
§ In the absence of the screening gate the plasmon speed s diverges at small wavevectors, see ((2]) at 
d — > oo. Although this effect suppresses the quantum fluctuations, it is not sufficient to restore the 
long-range order [28] , 
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nearest neighbor spins needs to be taken into account. The Hamiltonian describing the 
spin excitations takes the form 

F ff = ^JS r S Wl (4) 

i 

where Si is the spin at site I. As the exchange processes involve tunneling through a 
high barrier, the exchange constant is exponentially suppressed [34J, [35], 36J, 

J oc exp ( ) , (5) 



y/na B 

where r\ rs 2.80 [SB ED see also Sec. 14.1.11 Taken together, equations ([3]) and PJ 
account for all low-energy excitations of the one- dimensional Wigner crystal, i.e., the 
Hamiltonian of the system can be represented as the sum 

H = H p + H a . (6) 

Because of the absence of long-range order, one expects that in the low-energy limit the 
Wigner crystal should be a special case of the Luttinger liquid. The latter is commonly 
described [7J by a Hamiltonian of the form ([6]), with the charge and spin Hamiltonians, 
H p and H a , given by 

H P = J ^ [n 2 K p Ul + K' p \d x <P P f] dx, (7) 

Ha = I ^ [^ 2 K a Ul + K-\dM 2 ] dx + [ cos \V8Mx)] dx. (8) 

Here the bosonic fields p cr and n p (T describe the charge (p) and spin (a) excitations 
propagating with velocities u P)(7 . They obey canonical commutation relations 
[4> a (x),H a i(y)] = i5 aa /S(x — y). In the case of repulsive interactions, the Luttinger 
liquid parameter K p is in the range < K p < 1. The cosine term in (jBJ) is marginally 
irrelevant, i.e., the coupling constant gi± scales to zero logarithmically at low energies. 
At the same time, the parameter K a approaches unity as K a = 1 + gi±/27ru a . 

Both Hamiltonians ([3]) and (j7j) describe propagation of elastic waves in the medium. 
Their formal equivalence is established [36] by identifying 

7m y/2 T ^ irhn 

= U p( x ) = — rPl^)) u p = s, K p = - . (9) 

sj2 Ttnh 2ms 

On the other hand, even though both Hamiltonians @ and (|S]) describe spin excitations 
in the system, their equivalence is not obvious. Indeed, Hamiltonian (T4]) is expressed 
in terms of spin operators Si of the electrons, whereas its Luttinger-liquid analog (jSJ) 
is expressed in terms of the bosonic fields (f) a and U a . The connection is established 
via the well-known procedure [7] of bosonization of the Heisenberg spin chain ([41 . This 
procedure is applicable at energies much smaller than the exchange constant J, and 
reduces Hamiltonian (T4]) to the form (jSJ), see Ref. [36]. One therefore concludes that at 
low energies the Wigner crystal can indeed be viewed as a Luttinger liquid. 

It is important to point out, however, that the equivalence of the Wigner crystal 
and Luttinger liquid holds only at very low energies, e < J. Given the exponential 
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dependence (J5J) of the exchange constant on density, one can easily achieve a regime 
when an important energy scale, such as the temperature, is larger than J. In this case 
the bosonization procedure leading to (jHJ) is inapplicable, and the form (TJJ should be 
used instead. On the other hand, as long as temperature and other relevant energy scales 
are smaller than the Fermi energy, the charge excitations are bosonic and adequately 
described by either Hamiltonian (jHJ) or (J7J). 

2.2. Spin-charge separation in the one- dimensional Wigner crystal 

The Hamiltonian (JHJ)- (JHJ) of the Luttinger liquid consists of two separate commuting 
contributions associated with the charge and spin degrees of freedom. Consequently, 
the low-energy excitations of the system are charge and spin waves, decoupled from 
each other, and propagating at different velocities u p and u a . The operator annihilating 
a (right-moving) electron with spin 7 in this theory has the form 

^ Rl ^ = ^2twc ~ exp { ± ^7f^ CT ^ ~~ > ( 10 ) 

in which the charge and spin contributions explicitly factorize. (Here a is a short- 
distance cutoff, kp = im/2 is the Fermi wavevector of the electrons, and the +/— sign 
corresponds to electron spin 7 =t, !•) 

The Hamiltonian of the Wigner crystal (EJ) also consists of two commuting 
contributions describing the charge and spin degrees of freedom, with the main difference 
being the different form (j4]) of H a . However, the analogy with the Luttinger liquid is 
not complete, as the electron annihilation operator no longer factorizes [371 138] , 

Here the operator Zi n acts upon any state of the spin chain and produces a state 
with one less spin by removing spin 7 at site number /. The form of the fermion operator 
( TTTj) reflects the fact that when an electron is removed from the Wigner crystal, one of 
the sites of the spin chain (j3J) is also removed. In the absence of plasmon excitations, the 
sites are equidistant, and the site at point x has the number / = nx. On the other hand, 
if plasmons propagate though the crystal, the electrons shift by a distance proportional 
to (j) p , and the spin is removed from the site I = nx+^- (j> p (x), see f lTTj) . Thus the absence 
of factorization of the charge and spin components of the fermion operator (iTTj) reflects 
the simple fact that the spins Si in the spin chain (pQ) are attached to the electrons. 

The absence of spin-charge separation in the Hamiltonian of the Wigner crystal 
manifests itself if the system is not uniform, such as in the case of a quantum wire with 
a low electron density that depends on position, n = n(x). Assuming that the variations 
of n(x) occur at a length scale much larger than the distance between electrons, one 
can still bosonize the charge modes near every point in space, while accounting for the 
x-dependence of the parameters u p and K p . Thus one obtains 

H p = f {n 2 K p (x)Il 2 p + [K p {x)]-\d x <j> P ) 2 } dx. (12) 



^Rii x ) = 7= — ex P \ -/=[20p(aO - 9 p {x)\ \ Zi a (11) 

- 1 ■ l=nx+^4> p {x) 
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Figure 2. A quantum wire formed by applying negative voltage to the gates placed on 
top of a two-dimensional electron system. Electrons in the narrow channel between the 
gates are one-dimensional and their density is sufficiently low to achieve the Wigner 
crystal regime. Away from the center of the wire the electron density increases, and 
even the short range ordering of electrons is destroyed by quantum fluctuations. 



The exchange constant J in the Hamiltonian of the spin chain also acquires an 
x-dependence, as it clearly depends on the electron density, see flH]). Thus the spin 
Hamiltonian takes the form 



where %\ is the initial position of the Z-th electron. The appearance of the charge field 
4>p in H a again accounts for the fact that the plasmons shift the site I of the spin 
chain from its initial position by — p . Therefore the two contributions H p and H a to 
the Hamiltonian of the Wigner crystal commute only in the uniform system, when the 
exchange constant J does not depend on position. 

2.3. Conductance of a Wigner- crystal wire 

In experiment, the quantum wires are usually made by confining a two-dimensional 
electron system to a one-dimensional channel. One of the most common techniques is 
to place two metal electrodes above a GaAs heterostructure in which a two-dimensional 
electron system is formed, Figure [2J When a negative voltage is applied to the gates, 
the resulting electrostatic potential repels the electrons from the regions covered by 
the gates, but a narrow channel of electrons between the gates may still remain. The 
resulting quantum wire connects two large regions of two-dimensional electrons, which 
play the role of contacts to the wire. If the gate voltage V g is properly tuned, the electron 
density in the center of the wire can be sufficiently low for a Wigner crystal to form. 
On the other hand, the gates do not affect the electron density and the nature of the 
electron liquid in the two-dimensional leads. 

The physics of interacting electrons in two or three dimensions is very different from 
that of one-dimensional systems. Although at extremely low densities the electrons will 




(13) 
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form a Wigner crystal, this does not happen in typical GaAs heterostructures. Instead, 
the electrons are believed to be in a conventional Fermi-liquid state with quasiparticle 
excitations obeying Fermi statistics and carrying the charge of a single electron. In a 
one-dimensional system such a situation may only occur in the absence of interactions, 
as otherwise a Luttinger liquid state with bosonic excitations is formed. In the absence 
of interactions, however, the Fermi-liquid and Luttinger-liquid pictures are equivalent. 
Thus it is convenient to model the quantum wire device by a one-dimensional model 
with position-dependent interactions and electron density. In the central part of the 
system the density is small so that the interactions may be effectively strong. This 
region models the quantum wire. As one moves away from the central region, the 
density grows, the interactions become small, and asymptotically at large distances the 
electrons become noninteracting. These two semi-infinite noninteracting regions model 
the two-dimensional leads. 

Such a model was used in [9j [TOj [11] to calculate the conductance of a quantum 
wire described by the Luttinger liquid model. The Hamiltonian studied was essentially 
identical to f|T2|) . as the electrons were assumed to be spinless and only charge modes 
needed to be accounted for. It was demonstrated that the dc conductance of the wire 
is not affected by the interactions and remains quantized at e 2 /h. Let us illustrate this 
result with a simple semiclassical calculation. 

We start with the homogeneous wire, and for simplicity, instead of the Hamiltonian 
(J7|) we will use the equivalent form ([3]). Unlike papers |9j [101 [11], where a term was 
added to the Hamiltonian in order to describe the bias voltage applied to the wire at 
point x = 0, we consider a setup in which the wire is connected to a current source. A 
small ac current with frequency uj can be represented in terms of the velocity u of the 
elastic medium and the electron density as 

neii\ x= o = Iq coscut. (14) 

This expression should be viewed as a time-dependent boundary condition imposed on 
the elastic medium. As a result the medium begins to move periodically with frequency 
u, and plasmons propagating into the infinite leads dissipate power W = IqR/2 from 
the current source, where R is the resistance of the system. Let us calculate W in terms 
of the parameters of the elastic medium. Since the plasmons carry the energy of the 
oscillating medium in two directions at speed s, we can express the dissipated power as 

W = 2s(8), (15) 

where {£) is the energy density of the system. The latter consists of two contributions, 
the kinetic and potential energies represented by the two terms in (J3j). In a harmonic 
system the time-averaged values of the kinetic and potential energies are equal, so we 
will evaluate {£) by doubling the kinetic energy, 

(£) = mnii 2 = ^- J 2 (cos 2 ^) = -^-J 2 , (16) 

where we expressed the velocity u in terms of the current using (|T4|) . Substituting this 
expression into (fT5|) and comparing the result with the Joule heat law W = IqR/2, we 
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find the resistance 

R= 2jns = h_ ± 
e z n e z Vp 

where we used the density n = kp/ir for spinless electrons and defined the Fermi velocity 
in the interacting system as Vp = hkp/m. 

In the noninteracting limit, where the Luttinger liquid theory reproduces the low- 
energy properties of the Fermi gas, the plasmon velocity s = vp, and we recover the 
well-known result R = h/e 2 . The model considered in [9l [TUl [11] was described by 
the Hamiltonian ( 1T21) of the inhomogeneous Luttinger liquid, where the interactions are 
present only in a region of finite size L, modeling the wire, and vanish at x — > ±00. It 
is easy to see that the above calculation of the resistance is applicable to such a system 
as long as the low-frequency limit is considered. Indeed, at uo — > the wavelength of the 
plasmons ~ s/u is much larger than L, so the emission of the plasmons occurs in the 
noninteracting leads. Thus we have recovered the result [HI [El EH] for the conductance, 
G = e 2 /h. 

Our simple calculation also enables us to interpret the absence of corrections to the 
conductance due to electron-electron interactions in a finite region of a one- dimensional 
system. In the Luttinger liquid theory the main effect of the interactions is to change 
the compressibility of the electron system, thereby affecting the second term in ([3]). In 
the dc limit the wavelength of the plasmons is infinitely large, and thus the deformation 
d x u within the finite-size interacting region is negligible. Thus the system behaves as a 
noninteracting one. 

The above result for the spinless Luttinger liquid can be easily generalized to the 
case of electrons with spin. As we discussed in Sec. \2.2\ within the Luttinger-liquid 
approximation the charge and spin degrees of freedom are not coupled. Thus the applied 
bias or electric current couples only to the charge modes, and the above discussion can 
be repeated with the only modification being the different relation n = 2kp/n between 
the density and the Fermi wavevector. Substituting this expression instead of n = kp/rc 
in (fTTj) we find the resistance of the charge modes 

R P = ^ (is) 

and thus the expected doubling of the conductance, G = 2e 2 /h. 

On the other hand, we saw in Sec. 12.21 that in the inhomogeneous Wigner crystal 
there is no spin-charge separation, i.e., the Hamiltonian (fT3l) of the spin excitations 
depends explicitly on the charge field (j) p . One can therefore expect that the spin degrees 
of freedom will affect conductance when the Wigner crystal is not equivalent to the 
Luttinger liquid. Indeed, we show below that the spins have a significant effect on the 
electronic transport at temperatures T > J. 

In treating a one-dimensional Wigner crystal attached to noninteracting leads one 
has to overcome a fundamental problem caused by the lack of quantitative theory for 
the crossover regions that connect them, Figure [2j In the case of spinless electrons 
both the Wigner crystal and the leads can be viewed as special cases of the Luttinger 
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liquid, assuming that one is only interested in the low-energy properties of the system. 
Thus one can use the model ffl2|) of the inhomogeneous Luttinger liquid and obtain 
reliable results, provided that the exact form of the x-dependences of the parameters 
is not important. In the presence of spins there is an additional complication caused 
by the fact that the spin sector of a Wigner crystal is described by the Hamiltonian of 
a Heisenberg spin chain (jlj) because the spins are attached to well-localized electrons. 
Such a description is appropriate in neither the crossover region nor the leads, where 
the short-range crystalline order is absent. In our further discussion we will nevertheless 
use the model of the inhomogeneous spin chain (fl3l) for the whole system. This model 
is justified if the temperature is small compared to the Fermi energy in the center of the 
wire. When one moves away from the center, the density n grows, and consequently the 
exchange constant J rapidly grows, see (Q. Even if in the center of the system we had 
J<T, the crossover regime J ~ T will occur while the wire is still in the Wigner crystal 
regime, as J is still small compared to Ep. Eventually, when one moves sufficiently far 
from the center of the wire the exchange J becomes of order Ep, and the spin chain 
model is no longer appropriate. However, since in those regions we have J > T, the 
Heisenberg model (TjJ is equivalent to the spin sector ([8]) of the Luttinger liquid theory 
appropriate for both the crossover regions and the leads. Thus, at T Ep, one can 
describe the spin properties of the system by the model (fT3l) of an inhomogeneous spin 
chain as long as the exact shape of the dependence J (I) does not affect the results. 

Formally the quantum wire will be described by the Hamiltonian H p + H a given by 
( fl2l) and ( fl3l) . The electron density has a minimum at the center of the wire, resulting 
in an exponentially small exchange constant J, Figure El Far from the center of the wire 
the exchange constant reaches the value ~ E F . Since J depends on position, the spin 
excitations are coupled to the charge excitations. To find the resulting correction to the 
conductance of the wire, it is convenient to consider the setup of fixed current through 
the wire. Given the standard bosonization relation between d x <p p and the electron 
density, by fixing the current / at point x = one imposes the boundary condition 
4> p (0,t) = — [tx / \/2)q{t) on the charge modes, where q(t) is the charge transferred 
through the wire, i.e., / = eq. As discussed above, at small frequencies u the plasmon 
wavelength is very large, and electrons move in phase over distances much longer that 
the length of the wire. One can therefore replace <fi p (x) by its value at x = everywhere 
within the range where J depends on position, and convert the Hamiltonian (Fl3|) to the 
form 



The advantage of this form of H a is that it now commutes with H p . This does not mean 
that spin-charge separation is restored, as the spin excitations are still affected by the 
electric current. 

An immediate consequence of the commutativity of H p and H a is that the 
application of electric current through the wire gives rise to independent excitation 
of the charge and spin modes. Assuming that the power dissipated in each channel is 
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Figure 3. (a) The electron density as a function of position has a minimum in the 
center of the wire (x = 0), where nas <§C 1 and the Wigner crystal is formed. In the 
lead regions, nas is assumed to be large such that the interactions can be neglected, (b) 
The low density in the wire results in the exponential suppression ([5]) of the exchange 
constant J. In the lead regions J(l) saturates at Jqo ~ Ep. 



quadratic in current, we conclude W = IqR/2 = Iq(R p + Ro)/2. Thus the resistance of 
the wire is a sum, 



of two independent contributions due to the charge and spin excitations. Since we 
have already discussed the contribution (fl8|) of the charge excitations, we now turn our 
attention to R a . 

The spin contribution to the resistance depends crucially on whether the 
temperature is small or large compared to the value J of the exchange constant in the 
center of the wire, see Figure El At T <C J one can bosonize the spin excitations, i.e., 
convert H a to the form ([8]) with position-dependent parameters. Within this approach, 
an attempt to account for the coupling to the charge modes in ([13]) would result in 
corrections cubic in the bosonic fields. Such corrections are irrelevant perturbations, 
which are usually neglected as their contribution vanishes at T —>■ 0. Thus one concludes 
that R a — in the limit T/J — > 0. 

The absence of dissipation in the spin channel at low temperature can be interpreted 
as follows. The low-energy excitations of a Heisenberg spin chain are the so-called 
spinons [39J with spectrum 




where the wavevector k is defined in the interval (0,7r). At low temperature the state 
of the spin chain can be viewed as a dilute gas of spinons. Let us consider propagation 
of spinons in the spin chain ( Tl9|) with non- uniform J, Figure EJ^b), assuming for the 
moment q(t) = 0. If the variation of </(/) is very gradual, one can use the spectrum (I2T!) 
with Z-dependent exchange J. As a spinon propagates through the wire, its energy is 
conserved, but its momentum and velocity change because of the variation of J along 
the system. Clearly, if the energy of a spinon is less than ttJ/2, where J is the smallest 
value of the exchange constant in the system, Figure E^b), it passes through the wire 
without scattering. Conversely, spinons with energies exceeding irJ/2 are backscattered, 
Figure HI 



R — Rp + R, 



(20) 




e{k) = — sin k, 



(21) 
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Figure 4. Scattering of spinons at the quantum wire. Spinons with energies below 
7rJ/2 (shown in blue) slow down in the wire, but continue to move forward to the 
opposite lead. Spinons with energies above n J/2 (shown in red) stop before they reach 
the center of the wire and are scattered back. 



At q(t) 7^ the dependence J (I) shown in Figure MJo) is not static, but rather 
oscillates in position with respect to the spin chain. (More physically, the ac current 
moves the Wigner crystal with respect to the quantum wire, causing the time dependence 
of the exchange constants in (fl9i) .) The spinons passing through the wire without 
scattering are not affected by this oscillation. On the other hand, the spinons with 
energies e > nJ/2 are reflected by a moving scatterer. Such processes do change the 
energy of the spinons, and eventually lead to dissipation. At low temperature T< J 
the density of such (thermally-activated) spinons is very low, and one expects only an 
exponentially small resistance in this regime, 



It is worth mentioning that the resistance (l22i) is caused by excitations with energies of 
the order of the spinon bandwidth J. Such a correction cannot in principle be obtained 
by the bosonization procedure, which is accurate only at energies much smaller than J. 

The expression ( 1221) implies that the resistance R a grows with temperature. At 
T ^> J one expects this growth to saturate. Indeed, in this limit one can assume that 
J = in the center of the wire, i.e., the propagation of spin excitations through the wire 
is no longer possible. On the other hand, in the leads one still has T <C ~ Ep, and 
the picture of a dilute spinon gas still applies. Every spinon moving toward the wire is 
reflected back, resulting in a finite dissipation that no longer depends on J. 

Unfortunately, one cannot easily develop the theory of scattering of spinons in 
this regime, as such processes occur in the region where J (I) ~ T, and the spinon 
gas is no longer dilute. One can, however, conjecture that the dissipation resulting 
from all the spin excitations being reflected by the wire is universal in the sense that 
it does not depend on the exact nature of the scatterer. Thus if one can solve another 
problem where all the spin excitations in a one-dimensional system are reflected by a 
moving scatterer, the result for R a should be the same. The simplest example of such 
a problem is obtained in the same Wigner- crystal setup in the presence of a magnetic 
field B sufficient to polarize electrons in the center of the wire, T, J <C \ibB <C Ep, 
where //g is the Bohr magneton. Then only the electrons with spin directed along the 
field propagate through the wire whereas the electrons with opposite spin are confined 




(22) 
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to the leads. This problem can be easily solved in the framework of the bosonization 
approach [36], resulting in 

R* = T > J. (23) 

The result is easily understood by noticing that in combination with (1201) and (1T81) one 
finds the conductance G = e 2 //i which is the expected result for the conductance of a 
spin-polarized wire, where only one type of charge carriers participates in conduction. 
By our conjecture, the same reduction of conductance from 2e 2 /h to e 2 /h occurs in the 
absence of the field, provided J<T, because in both cases all the spin excitations are 
reflected by the wire, resulting in the same dissipation. This conclusion is consistent 
with some of the measurements of the conductance of quantum wires a low density 
[HI El El E], showing a small plateau at G = e 2 /h. 

3. Classical transition to the zigzag structure 

In section [21 we discussed the physics of a purely one-dimensional crystal. 
Experimentally, however, quantum wires are created by confining three-dimensional 
electrons to a narrow channel by an external confining potential. The electron system 
in the wire can be viewed as one-dimensional as long as the typical energy of the 
transverse motion is large compared with all other important energy scales; otherwise, 
deviations from one-dimensionality arise. The remainder of this review addresses 
the resulting quasi-one- dimensional physics, starting with the classical transition from 
a one-dimensional to a quasi-one-dimensional Wigner crystal that was studied in 
Refs. [301 EH [32]. 

To be specific, we consider here a confining potential that mimics the experimental 
situation. In a typical setup the confining potential in one direction, say the ^-direction, 
is provided by the band bending at the interface of two semiconductors with different 
band structure (typically GaAs and AlGaAs). This provides a very tight confinement 
and, correspondingly, the energy scales for transverse excitations are large. Therefore, 
at low energies, the possibility of electron motion in the ^-direction may be neglected. 
By contrast, confinement in the ^-direction is provided by nearby metallic gates which 
create a relatively shallow confining potential. Deviations from one- dimensionality arise 
due to lateral displacements in this shallow potential which may be assumed parabolic: 

Konf = \mn 2 Y,vl (24) 

i 

where Q is the frequency of harmonic oscillations in the confining potential, and yi is 
the transverse coordinate of the electron at site i. 

As the electron density n grows, so does the typical energy VJ nt ~ (e 2 /e)n of the 
Coulomb interaction between electrons. Eventually, it becomes energetically favorable 
for electrons to move away from the axis of the wire. This happens when the distance 
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between particles is of the order of the length scale 



r 



2e 2 



(25) 



emQ 2 ' 

defined by the condition that the confinement and the Coulomb repulsion, K on f(r ) = 
|mf2 2 rQ and Vi nt (r ) = e 2 /er , are equal [32] . 

The quasi-one-dimensional arrangement that maximizes the distance between 
electrons — and consequently minimizes the Coulomb interaction energy VJ nt = 
(e 2 /e) J2i<j \ r i~ r j\^ 1 — a t a given cost of confining potential energy is a zigzag structure, 
see Figure [T^b). The exact shape of the zigzag crystal can be found by minimizing its 
energy per particle 
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(26) 



with respect to the distance w between the two rows of the zigzag crystal. Here v = nr 
is the dimensionless density, the first two terms account for the interactions between 
electrons within the same row and in different rows of the zigzag structure, respectively, 
and the last term stems from the confining potential. 

One finds that the distance between rows is given by the solution of the equation 



v 

T 



1 



V 



- 1 [(i-i) 2 + ^ 

Below the critical density [30j EI] 



3/2 



W 



0. 



(27) 



7C(3) 



0.780, 



(28) 



the only solution is w = and, therefore, the crystal is one- dimensional. At densities, 
v > u c , a lower-energy solution with w ^ appears, and the zigzag structure is 
formed. The distance between the two rows of the zigzag crystal grows with density. 
In particular, just above the transition point u c , the distance between rows behaves as 
ro[-\/24/93C(5)/^ 2 ]A/5i/, where bv = v—v c . Upon further increasing the density, the 



w 

zigzag crystal eventually becomes unstable at v ~ 1.75. At larger densities, v > 1.75, 
structures with more than two rows are energetically favorable [32] . 

Such a classical description of the system is valid only in the limit where the distance 
between electrons is much larger than the Bohr's radius, n^ 1 ^> ag. As the zigzag regime 
corresponds to distances between electrons of order r , it can only be achieved if r is 
sufficiently large, tq ^> as- This motivates the introduction of a density- independent 
parameter 

(29) 



rn = — , 
a B 



which characterizes the strength of Coulomb interactions with respect to the confining 
potential. If tq <C 1, as the electron density grows, the interactions become weak at 
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n ~ a^ 1 C r _1 . As a result, the one-dimensional Wigner crystal melts by quantum 
fluctuations before the zigzag regime is reached. By contrast, if ^> 1, interactions 
are still strong [nas <C 1) at densities n ~ r$ , and the classical description of the 
transition to the zigzag regime is applicable. As oc fi~ 2 / 3 , the strongly interacting 
case therefore requires a shallow confining potential. Note that the condition 3> 1 
can be rewritten as W 3> Ob, where W = ^Jh/mVl is the (quantum) width of the wire. 



4. Spin properties of zigzag Wigner crystals 

In a Wigner crystal electrons are localized near their lattice positions due to the mutual 
Coulomb repulsion. The potential landscape thus created is such that any deviation from 
these lattice positions incurs an increase in Coulomb energy. In particular, the exchange 
processes which give rise to spin-spin interactions require tunneling of electrons through 
the Coulomb barrier that separates them. As pointed out in section 12.11 the resulting 
spin couplings in a one-dimensional crystal are fairly simple: as the tunneling amplitude 
decays exponentially with distance, only nearest neighbor exchange processes have to 
be taken into account. Thus, the spin degrees of freedom of a one-dimensional Wigner 
crystal are described by an antiferromagnetic Heisenberg chain (j3j) with nearest neighbor 
exchange energy J whose properties were discussed in section |2j 

In a zigzag chain, spin couplings become more interesting. Close to the zigzag 
transition, the nearest neighbor exchange is dominant as in the one-dimensional case. 
However, as the zigzag structure becomes more pronounced each electron is surrounded 
by four close neighbors rather than only two as in the one- dimensional crystal, and, 
therefore, the next-nearest neighbor couplings can no longer be neglected. Instead of 
one coupling constant, one needs to take into account a nearest neighbor exchange 
constant J\ and a next-nearest neighbor exchange constant J<i- Both couplings are 
antiferromagnetic and, therefore, compete with each other. If J 2 is large enough 
(J2 > 0.24.. Ji [HH HH H2]), the antiferromagnetic ground state gives way to a dimer 
phase characterized by a non-vanishing order parameter D oc ((S^i+i — -S^i-i) • S 2 i) 
and a resulting spin gap. The dimer structure is particularly simple on the so-called 
Majumdar- Ghosh [4l3l \M\ line J 2 = 0.5Ji, where the dimers are just nearest neighbor 
singlets. The magnitude of the spin gap is a non-monotonic function of the ratio Jij J\. 
it reaches its maximum close to the Majumdar- Ghosh line and becomes exponentially 
small at J 2 J\. 

It turns out, however, that these two-particle exchanges are not sufficient to describe 
the spin physics of the zigzag crystal. In addition, ring exchanges, i.e., cyclic exchanges 
of n > 3 particles, have to be taken into account. Defining exchange constants in such 
a way that they are all positive, the Hamiltonian of the system then reads 

-^ring = 2 (y^W + ^2-P/H-2 — -hi.Pu+ll+2 + Pl+2l+ll) 

I 

+ Ji{Pimwv& + Pmi-\zi+ii) — • • • ) , (30) 
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where P^ is a permutation operator and Pi x ...i N = Pixi 2 Pi 2 i 3 ■ ■ ■ P% N ix- Here we still label 
particles according to their position along the wire axis only: thus, nearest neighbors 
are particles in opposite rows whereas next-nearest neighbors are the closest particles 
within the same row. Note that for densities in the range 1.45 < v < 1.75 the lateral 
displacement w is so large that the distance between nearest neighbors becomes larger 
than the distance between next-nearest neighbors. 

Ring exchanges are interesting because they might stabilize a ferromagnetic ground 
state. While exchanges involving even numbers of particles favor a spin-zero ground 
state, exchanges involving odd numbers of particles favor a ferromagnetic arrangement 
of spins [45]. Thus, the simplest ring exchange process that could lead to a polarized 
ground state is the three-particle exchange. In fact, ring exchanges have been extensively 
studied in two-dimensional Wigner crystals [46j H3, HU H9]. In that case the three- 
particle ring exchange dominates in the low-density limit which implies a ferromagnetic 
ground state of the strongly interacting Wigner crystal in two dimensions]]]] To find out 
whether the physics of the zigzag Wigner crystal is similar, one needs to compute the 
exchange constants for nearest neighbor, next-nearest neighbor, and the various ring 
exchanges. 



4-1. Computation of exchange constants 

To introduce the method, we start by discussing the one-dimensional case where the 
only non-negligible exchange is the nearest neighbor exchange. 



4-1.1. Exchange constants for the one- dimensional Wigner crystal The nearest 
neighbor exchange constant J can be determined by computing the tunneling probability 
of two electrons through the Coulomb barrier that separates them. If the barrier 
is sufficiently high and, therefore, tunneling is weak, one may use the semiclassical 
instanton approximation. This corresponds to finding the classical exchange path in the 
inverted potential by minimizing the imaginary-time action. 

It is convenient to rewrite the action in dimensionless form by rescaling length in 
units of 1/n and time in units of a/ em/e 2 n 3 . The action of the system is then given as 



h f 

Sid = . Vid, where Vm[{xj(r)}] = / dr 

x/na B J 



x 2 ^ 1 



(31) 

j j<i 

As a first approximation one may fix the positions of all particles except the two that 
participate in the exchange process, say j = 1 and j = 2. Symmetry fixes the center of 
mass coordinate of the exchanging electrons and, therefore, the minimization has to be 
done only with respect to the relative coordinate x = x<i — X\. The tunneling lifts the 
ground state degeneracy present due to inversion symmetry x — > —x, and the exchange 
energy can be identified with the resulting level splitting. 

|| The ferromagnetic state is predicted to occur only at extremely low densities characterized by a value 
of r s > 175 |49j . where r s is the ratio of the Coulomb interaction energy to the Fermi energy. 
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(a) (b) 




Figure 5. Sketch of typical exchange paths for (a) v -C v c and (b) v < v c . the size 
of the loop where electrons move away from the axis of the wire is determined by the 
length scale tq. 



The instanton approximation yields the exchange constant J in the form ([5]), where 
7] is the dimensionless classical action obtained from the minimization procedure. One 
finds r] « 2.817 [50]. At low densities, nas <C 1, the exponent is large leading to 
exponential suppression of J, and thus the prefactor omitted in is of secondary 
importance. 

Fixing the positions of all particles except the two participating in the exchange 
process is a somewhat crude approximation. Neighboring electrons see a modified 
potential due the motion of the exchanging particles and, therefore, experience a force 
that displaces them from their equilibrium positions. A better estimate for r) can 
be obtained by including these mobile "spectator" particles in the minimization. By 
allowing spectators to move during the exchange process, one expects to find a reduced 
value for rj because more variables are varied in the minimization procedure. It turns 
out, however, that the effect is very small. As more spectators are added, rj approaches 
the asymptotic value rj ~ 2.798 [50], EU [52] , i.e., the result changes by less than 1%. 

4-1.2. Exchange constants for the zigzag Wigner crystal In the presence of a confining 
potential, the motion of the exchanging electrons is no longer restricted to one dimension, 
i.e., the position of an electron is now given by a two-dimensional vector rj = (xj,yj). 
In particular, if the wire width W is larger than the Bohr's radius a# or, equivalently, 
the interaction parameter introduced in Eq. ( 1291) is large, rn ^> 1, electrons can make 
use of the transverse direction to go "around" rather than "through" each other during 
the exchange process. This reduces the Coulomb barrier and, therefore, increases the 
tunneling probability. The characteristic length scale of the transverse displacement 
is given by the length r , introduced in section [3] Typical trajectories for the one- 
dimensional crystal are shown in Figure [5] for low and moderate densities, v <C v c and 
v < u c , respectively. At low densities, v <C u c , the exchange part follows the bottom of 
the confining potential until electrons come within a distance of order tq of each other. 
Thus, only a small part of the exchange path explores the transverse direction, leading 
to a relatively small correction to the tunneling action S\jy. The results of section 14.1.11 
are recovered in the limit v — > 0. As one approaches the transition to the zigzag crystal, 
the exchange trajectories become more and more two-dimensional and consequently the 
exchange couplings are modified significantly. Finally, at v > u c , also the equilibrium 
positions of the particles are displaced in the ^/-direction. 

The exchange constants for the zigzag Wigner crystal can be obtained in the same 
way as for the one-dimensional Wigner crystal [53]. However, by contrast to the one- 
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dimensional case, the structure of the zigzag crystal chang function of density. 

As a consequence the rescaling of lengths and times used in the one-dimensional case 
is not appropriate here. A dimensionless action in a transverse confining potential is 
conveniently defined using the interaction parameter r^. Namely 



S 2 d = ^V^V2D, where r?2D[{rj(r)}] = / dr 



3 J<* 



(32) 



Here lengths have been rescaled in units of ro whereas times has been rescaled in units of 
y/2/Q. Furthermore, comparing fl3T|) and fl32|) . the differences are the one-dimensional 
vs two-dimensional coordinates and the additional term due to the confining potential 
in (J32]). 

As a result the exchange constants take the form 

J i = J i exp (-TftV^fi) , (33) 

where J\ is the nearest-neighbor exchange constant, J 2 is the next-nearest neighbor 
exchange constant, and J\ for / > 3 is the exchange constant corresponding to the l- 
particle ring exchange. The exponents rji are obtained by minimizing the dimensionless 
action '72D[{ r j(' r )}] for a given exchange process. Whereas in the strictly one-dimensional 
case 77 was just a number, now the electron configuration changes as a function of density 
and, therefore, the exponents r]i depend on density, too. 

Note that while in the one- dimensional case the inclusion of spectators had little 
effect on the results, here the spectators turn out to be much more important [53] . Figure 
My^) shows the change of the exponents rji as spectators are included. As one can see, 
the first few spectators modify the results significantly. However, the results converge 
rapidly as more and more spectators are added. Thus, the spin couplings are generated 
by processes that involve the motion of a small number of close-by electrons. Therefore, 
these couplings should not be affected by deviations from the perfect crystalline order 
at large distances, caused by quantum fluctuations. 

Figure E](b) shows the exponents rji as a function of the dimensionless density 
v [53| 155] . Ring exchanges with more than four particles are not included as they are 
negligibly small at all densities. At small v < 1.2, the crystal geometry is still close to 
one-dimensional. In that regime r]i is the smallest exponent and therefore, as expected, 
the nearest-neighbor exchange J\ dominates. However, as density increases and the 
distances between nearest neighbors and next-nearest neighbors become comparable, in 
the regime 1.2 < v < 1.5 the three-particle ring exchange constant J3 becomes largest. 
Finally, at even higher densities v > 1.5 the four-particle ring exchange is dominant 
(until the zigzag crystal gives way to structures with more than two rows at v « 1.75). 
In the next section the ground states generated by these spin couplings will be discussed. 



4-2. Spin phases of the zigzag Wigner crystal 

In order to extract the spin properties of the ground state, it is convenient to rewrite 
Hamiltonian ( 1301) in terms of spin operators using the identity = | + 2 Si ■ Sk- In 
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Figure 6. (a) Dependence of the exponents i]i on the number of spectators included 
in the calculation [54]. Results are shown for v = 1.3. (b) Exponents r\i for the nearest 
neighbor, next-nearest neighbor, three-particle ring, and four-particle ring exchange as 
a function of dimensionless density v [55] . 



the absence of ring exchanges the system is described as a Heisenberg spin chain with 
nearest neighbor and next-nearest neighbor coupling, 

#12 = ( Jl Sl ■ + J2 Sl ■ Sl +^ ■ ( 34 ) 

I 

As discussed at the beginning of this section, depending on the ratio of J\ and J 2 , one 
finds an antiferromagnetic and a dimer phase. The contribution of the three-particle 
ring exchange is 

H 3 = -J 3 ^ (2«i • S l+1 + S t ■ S l+2 ) . (35) 
1 

Thus, no new terms are generated — the Hamiltonian retains the same form (13^1) . albeit 
with modified coupling constants 

J X = J X - 2J 3 , J 2 = J 2 - h- (36) 

The important consequence is that the new coupling constants J\ and J 2 may now 
be either positive or negative, corresponding to antiferromagnetic or ferromagnetic 
interactions, respectively. The phase diagram of a Heisenberg spin chain with both 
antiferromagnetic and ferromagnetic nearest and next-nearest neighbor couplings has 
been widely studied in the literature [101 HU S21 SH SH [561 EZl EH ESI [60l EI]. In 
addition to the antiferromagnetic and dimer phases existing for positive couplings, a 
ferromagnetic phase appears. The phase diagram is shown in Figure [7(a). 

This phase diagram is sufficient to determine the ground state of the strongly- 
interacting zigzag Wigner crystal at low and intermediate densities. At low densities, 
the system is in the antiferromagnetic phase (Ji > 0, | «y 2 1 ^ J\)- At intermediate 
densities, the three particle ring exchange dominates. As a result both coupling 
constants become negative, Ji, J 2 < 0, and therefore the system is in the ferromagnetic 
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Figure 7. (a) Phase diagram of the Heisenberg spin chain with nearest neighbor 
coupling Ji and next-nearest neighbor coupling Ji- (b) Preliminary phase diagram 
of the zigzag spin chain including four-particle ring exchange J4, obtained by exact 
numerical diagonalization of finite-size chains [55] . When J4 is large, novel phases 
appear. (Triangles, squares, and circles correspond to the boundaries obtained for 
iV=16, 20, and 24 sites, respectively.) 



phase. The spontaneous spin polarization suggested as a possible explanation of the 
0.7 anomaly can, thus, occur in strongly interacting quantum wires, if deviations from 
one-dimensionality are taken into account. 

At higher densities, the situation becomes more complicated. While the three- 
particle ring exchange only modifies the nearest neighbor and next-nearest neighbor 
exchange constants, the four-particle ring exchange generates new terms in the 
Hamiltonian, namely a next-next-nearest neighbor exchange and, more importantly, 
four-spin couplings. The corresponding spin Hamiltonian reads 

l n=l 

+ 2 [(Si ■ Si + i)(Si + 2 ■ S1+3) + (Si ■ Si + 2)(Si + i ■ S1+3) — (Si ■ Si + s)(Si + i ■ Si+z)] J . 

The phase diagram in the presence of these couplings is not yet fully understood. First 
results were obtained using exact diagonalization of short chains with up to iV = 24 
spins [55]. If J 4 <C \ Ji\, I J 2 \, the same phases as in the Heisenberg spin chain without 
four-particle ring exchange appear as can be seen in Figure^b). However, as J4 becomes 
of the same order as the other coupling constants new phases appear. The simplest one 
to identify is a partially polarized phase (labeled £ M' in Figure 0(b)) adjacent to the 
ferromagnetic phase. While this phase seems to persist in size and shape as the number 
of spins increases, it is currently unclear whether it survives in the thermodynamic limit. 
In addition a region (labeled '4P' in Figure [3(b)) where the ground state is unpolarized 
but different from the antiferromagnetic and dimer phases occurs. This '4P' region could 
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Figure 8. Spin ground states of interacting electrons in quantum wires in the zigzag 
regime |55j . 



correspond to a single or several phases. Unfortunately, the size dependence in this part 
of the phase diagram turns out to be very complicated. Due to frustration introduced 
by the four-particle exchange, a large number of low-energy states exist. Therefore, the 
study of short chains does not allow one to determine the properties of the ground state 
in this regime. 

4-3. Spin phases of interacting quantum wires in the quasi- one- dimensional regime 

While the above results were obtained in the limit m ^> 1, interaction parameters in 
realistic quantum wires vary widely, ranging from tq < 1 in cleaved-edge overgrowth 
wires to « 20 in p-type gate-defined wires [62, 63, 64|. While the former are weakly 
interacting, the latter are clearly in the strongly interacting regime. However, the above 
analysis based solely on exponents is not sufficient to determine the ground state of 
interacting electrons in a quantum wire at finite r^. In order to obtain a phase diagram 
in that case, the prefactors have to be computed which can be done by including 
Gaussian fluctuations around the classical exchange paths. 

Using the exchange constants Ji(v, tq) computed in this way [55] and the phase 
diagram shown in Figure the ground states realized for given system parameters can 
be determined. The resulting phase diagram is shown in Figure It turns out that the 
partially and fully polarized phases are realized only at large > 50. At moderately 
large the transition occurs directly from the antiferromagnetic phase to a phase 
dominated by the four-particle ring exchange. These findings, thus, do not support 



the interpretation [ISl HH HH HSl [13 HH1 [201 EH [221 [231 [2H [25] of the so-called 0.7 



anomaly in terms of spontaneous spin polarization. 

Even at sufficiently strong interactions, the question arises of how a ferromagnetic 
state in the quantum wire manifests itself in the conductance. It is tempting to assume 
that in the fully polarized state the wire supports only one excitation mode and thus 
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Figure 9. Simple model of exchange coupling in a ferromagnetic Wigner crystal 
coupled to non-magnetic leads. The coupling constant vanishes in the contact region 
at points —a and a, and therefore the spin excitations in the leads and in the wire 
decouple. 



has conductance e 2 /h. This is indeed the case when the full polarization of electron 
spins is achieved by applying a sufficiently strong magnetic field. Such a field creates 
a gap in the spectrum of spin excitations, and below the gap the system is equivalent 
to a spinless electron liquid with conductance e 2 /h. It is important to stress that the 
situation is very different if the full spin polarization is achieved due to internal exchange 
processes in the electron system, rather than the external magnetic field. In this case, the 
ground state is degenerate with respect to spin rotations, and thus the system supports 
gapless spin excitations — the magnons. As a result, the conventional argument in favor 
of conductance value e 2 /h no longer applies. 

In studying the conductance of a ferromagnetic wire it is important to keep in 
mind that the properties of the electron system inside the quantum wire in general do 
not fully determine its conductance. Indeed, since the electric current flows between 
non-magnetic leads through a ferromagnetic wire, the spatial non-uniformity of the 
system needs to be considered carefully, and the problem of determining the conductance 
complicates considerably. In the case of a ferromagnetic zigzag Wigner crystal in the 
middle of the wire, the weakening of the confining potential in the contact region would 
lead to either melting of the crystal or the emergence of a crystal with more and more 
rows. In both cases modeling of the spin interactions in the transition region is by no 
means obvious. 

The simplest model that might capture the relevant physics is one where the system 
is described by Hamiltonian (J4j) with an effective position-dependent nearest-neighbor 
exchange constant J(x) as depicted in Figure [9j In the leads, interactions are weak and 
antiferromagnetic and therefore J is large and positive. In the wire, interactions are 
strong and ferromagnetic and therefore J is small and negative. Through the contact 
regions, J varies smoothly and changes sign at points —a and a. Within this model, the 
arguments of section 12.31 lead to the conclusion that the spin polarization does suppress 
the conductance. Namely, since the exchange coupling constant vanishes at the borders 
of the ferromagnetic region, i.e., at ±a, the spin degrees of freedom in the leads are 
decoupled from those in the wire and, thus, the propagation of spin excitations through 
the wire is blocked. Accordingly, the value of the conductance is reduced by a factor 
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2. By contrast to the antiferromagnetic case in one- dimensional wires, this suppression 
would persist down to temperatures T — > due to the vanishing of the spin coupling in 
the contact region. 

Of course the contact region in real quantum wires is more complex and a 
satisfactory theory for the conductance of strongly interacting quasi-one-dimensional 
quantum wires that correctly takes into account the spin degrees of freedom is an open 
problem. 

5. Orbital properties of zigzag Wigner crystals 

In addition to the spin physics discussed in the previous section, quasi-one-dimensional 
wires have interesting orbital properties. In this section, we discuss the transition from 
a one-dimensional to a quasi-one-dimensional state for the case of spinless electrons, 
based mainly on Ref. [65] . 

As shown in section [3], the classical transition from a one- dimensional to a zigzag 
Wigner crystal can be obtained simply by minimizing the energy of the interacting 
electron system in a transverse confining potential. A different way of studying the same 
transition is by considering the phonon modes of the crystal. In the one-dimensional 
crystal one longitudinal and one transverse phonon mode exist. The longitudinal phonon 
is gapless because sliding of the entire crystal along the wire axis does not cost any 
energy. On the other hand, the transverse mode is gapped with a gap frequency equal 
to the frequency of the confining potential. The transition to a zigzag state is driven 
by a softening of the transverse phonon at wave vector k = irn which corresponds to a 
staggered displacement of electrons transverse to the wire axis. In the zigzag crystal, 
we obtain two longitudinal and two transverse phonon modes with the following low-g 
dispersions [66] close to the transition {bvjv <C 1): 

cj|| (g) = ^nJul\n^ \q\, u^q) = V2SI + 0(q 2 ), (38) 

u±0 (q) = n + 0(q 2 ), cu ±n (q) = V6qJ— + q\ (39) 

V v c 48 

where q = kj (jrn). 

Thus, only at the transition point Su = 0, two gapless phonon modes exist with 
dispersions u>\\o(q) and uj^{q) = (7r/2\/2) VL^Jvf\n2 \q\. Within the zigzag regime, there 
is a single gapless excitation corresponding to in-phase longitudinal motion of the two 
rows that constitute the zigzag crystal. The soft-mode uj± n (q) describing the out-of- 
phase transverse motion acquires a gap A c \ oc y5u. 

This behavior is markedly different from the noninteracting case. In a 
noninteracting system the transition from a one-dimensional to a quasi-one-dimensional 
state happens when the chemical potential is raised above the subband energy of the 
second subband of transverse quantization. In the quasi-one-dimensional state, the 
two occupied subbands are decoupled and each of them supports a gapless electronic 
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Figure 10. At zero and infinite interaction strength the quasi-one-dimensional system 
supports two and one gapless excitation mode, respectively. Here the interaction 
strength is characterized by the parameter rn introduced in Sec. S Possible phase 
diagrams consistent with these findings arc shown: (a) A tricritical point exists at 
a finite interaction strength, (b) Weak quantum fluctuations destroy the gap of the 
classical Wigner crystal, (c) Already infinitesimally weak interactions induce a gap in 
the second mode. 

excitation mode, i.e., above the transition two gapless modes exist rather than just one 
as in the classical Wigner crystal. 

One might, thus, expect that the phase diagram of the system as a function 
of interaction strength is as shown in Figure [TOT a). namely two distinct quasi-one- 
dimensional phases exist at weak and strong interactions. Consequently one should 
find a tricritical point at a finite interaction strength where the nature of the transition 
from a one- dimensional to a quasi-one-dimensional state changes. However, the phase 
diagram Figure [TOT a) is not the only one consistent with both the above findings for 
the noninteracting case and the classical Wigner crystal at infinite interaction strength. 
Two alternatives are shown in Figure fTOT b.c). To distinguish between the different 
possibilities, one needs to study the nature of the transition as a function of interaction 
strength. In particular, the following questions have to be answered: (i) Do weak 
quantum fluctuations destroy the gap found in the classical zigzag crystal at strong 
interactions? (ii) Do infinitesimally weak interactions lead to a gap in the second mode 
above the transition? 

5.1. Quantum theory of the zigzag transition 

Let us consider the strongly interacting case first and account for the quantum nature of 
the system. In particular, using the classical Wigner crystal configuration as a starting 
point, we now include quantum fluctuations. The phonon modes ( l38l) and (J39l) reflect 
the fact that there are only two types of possible low-energy excitations: the longitudinal 
plasmon mode and a staggered transverse mode. It turns out that in the vicinity of the 
transition the two modes decouple. The acoustic spectrum of the plasmon mode is 
protected by translational invariance and, thus, at least one gapless excitation mode 
exists in the system. More interesting is the staggered transverse mode. 




Figure 11. Mapping of the <^ 4 -theory to a spin chain. 



To describe the transverse displacements of the electrons, a staggered field (pi = 
( — l) l yi is introduced. In the vicinity of the transition, ipi is slowly varying on the scale 
of the inter-electron distances and, therefore, the continuum limit <pi — ► <p(x) can be 
taken. Expanding the action up to fourth order in (p, one finds 

S[(p] = Ah^ J dr dx [(d T ip) 2 + {d Xr of -8vip 2 + <p A ] , (40) 

where the variables have been rescaled such as to provide the simplest action possible. 
The form of the action as well as all following conclusions do not depend on the exact 
shape of the confining potential. For a parabolic confining potential the constant A is 
given as A = [7C(3)] 3 / 2 v / hT2/(31C(5)). 

The classical transition point as discussed above corresponds to 8v = 0. Here 
the transverse mode becomes unstable, and the quartic term is needed to stabilize the 
system. Quantum fluctuations may affect both the transition point and the nature of 
the transition. A convenient way to analyze the quantum-mechanical problem is to 
refermionize. As a first step, we rediscretize the coordinate x along the wire axis. The 
discrete version of the Hamiltonian then describes a set of particles moving in a double- 
well potential Vow ~ ~^ l p] + ( as depicted in Figure [IT]) and interacting through a 
nearest-neighbor interaction oc (ipj — pj + \) 2 . If the double-well potential is sufficiently 
deep, i.e., if A is sufficiently large, the particles are almost completely localized in one 
of the wells at tpj = ±a/ A/2. Then each particle can be described by a pseudo-spin 
operator, namely (pj = a/ A / 2 erj , where cr| is a Pauli matrix. In terms of these new 
variables, the Hamiltonian consists of two terms corresponding to tunneling between the 
two wells and the nearest-neighbor interaction, respectively. Tunneling is described by 
H t = —t J2j a j whereas the interaction term reads H NN = —v J2j a j cr j+i- The resulting 
Hamiltonian H t + Hnn is the Hamiltonian of the transverse field Ising model [67] . Here 
the parameters t and v are related to the parameters of the original model. In particular, 
they can be tuned by changing the chemical potential which controls the transition, i.e., 
t = t(ji) and v = v(fi). 

In order to arrive at a fermionic description, a Jordan- Wigner transformation is 
used. It turns out that the Hamiltonian takes a much simpler form if one rotates 
<j x — > —<j z first. The representation of the spin operators in terms of (spinless) fermions, 

af = a 1 : ■ }a] = 2a)e™^<3 °H aj = cr]-}o] = 2e~' m ^ a ^ aj , a] = 2a] aj -l, (41) 

where aj and aj are fermion creation and annihilation operators on site j, respectively, 
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Figure 12. After the Jordan- Wigner transformation, one obtains a tight binding 
model for spinless fermions with bandwidth 2v and chemical potential — It. The 
transition from a one-dimensional to a quasi-one-dimensional state happens when the 
chemical potential reaches the bottom of the band. 



then yields the noninteracting Hamiltonian [6j 

Hi = ^2 [2* a ] a j ~ v ( a j~ a j S j ( a j+i + a j+i) • (42) 



Note that this Hamiltonian is essentially the transfer matrix of the two-dimensional 
classical Ising model [68J near the transition. The connection can be made clear 
by considering the mapping between d- dimensional quantum and (d + l)-dimensional 
classical models [69]. Thus, the one-dimensional quantum Ising model studied here is 
equivalent to the two-dimensional classical Ising model |70j . 

In Hamiltonian (142p one can identify three different contributions: — v (ajaj +1 + 



a \+i a j) describes a tight-binding model with bandwidth 4v, the local term 2ta\dj yields 
' ,t„t ' 



the chemical potential fMf = —2t of the spinless fermions, and finally — w(a]-a] +1 — o 3 -o 3 - + i) 
is a BCS-like pairing term with p-wave symmetry. The one-dimensional regime of our 
original model corresponds to t > v when the chemical potential lies below the bottom of 
the tight-binding band and, therefore, all the fermionic states are empty. The transition 
to a quasi-one- dimensional regime happens when the chemical potential reaches the 
bottom of the band at t — v . For t < v , some of the fermionic states describing the 
motion of the original electrons transverse to the wire axis are filled. Due to the pairing 
term in the Hamiltonian, these states acquire a gap. With the help of a Bogoliubov 
transformation, the Hamiltonian can be diagonalized to obtain the energy spectrum. 
As a result one finds that the gap is given as A = 2\t — v\. 

In terms of the original system parameters, one expects that t and v are non-singular 
functions of the chemical potential. The critical chemical potential /i c is thus defined 
by the condition t(/i c ) = v(fi c ). Furthermore, the behavior of the gap is obtained by 
expanding t — v in the vicinity of the transition point \x c . As a consequence one obtains 
a linear gap, 

Aoc|/i-/i c |. (43) 
While quantum effects modify the nature of the transition, the number of gapless 
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excitations remains the same: in the strongly interacting system, only one gapless 
excitation exists. Thus, the phase diagram of Figure [TDT b) is ruled out — weak quantum 
fluctuations do not destroy the gap of the transverse mode in the quasi-one-dimensional 
regime. To differentiate between the phase diagrams of Figures ITDT a) and ITOlfc) . a 
complementary approach has to be used. In the next section we consider the limit of 
weak interactions to check whether or not they lead to the formation of a gap just above 
the transition, as suggested by the scenario of Figure ITOTc). 



5.2. Two-subband system at weak interactions 

The limit of weak interactions can be treated using a renormalization group approach. 
Here the description in terms of two subbands due to transverse size quantization is 
a good starting point. Each subband is described by fermionic operators ipj, where 
j = 1,2. The free Hamiltonian is just 

-^J^+e^J, (44) 



where Sj are the subband energies. For a parabolic confining potential, Sj = hfl(j — ~). 

Interactions can be separated into intra-subband and inter-subband interactions. 
One needs four interaction constants to describe the system: gj ~ Vjj(0) — Vjj(2k F j) 
which describes intra-subband forward scattering, 

g x ~ MO) - \(V£(k F1 - k F2 ) + V£(k F1 + k F2 )) (45) 

which describes inter-subband forward scattering, and 

g t ~ V£(k F1 - k F2 ) - V™(k F1 + k F2 ) (46) 

which describes transfer of two particles between the subbands as shown in Figure [TBT a). 
Here 

Va(k) = jdxdx'e ik ^ jdydy'V iQt (r-r')xUy)x%y'), (47) 

V^{k) = j dxdx'e ik ^ J dydy' V int (r-r') Xl (y)xAy)My')xAy')^ (48) 

where Xi are the transverse eigenmodes in the confining potential. 

It is well known that forward scattering in one dimension does not open a gap in 
the system [7]. By contrast, the two-particle transfer between subbands described by 
the coupling constant g t could open a gap. To assess whether this is indeed the case, 
the renormalization group (RG) is used, i.e., reducing the bandwidth from D down 
to D, the scale-dependent coupling constants are determined. The RG equations for a 
two-band system are given as [TlJ [72j [73] 

1 1 1 

9i = - 9t, 92 = -7T^ 9t, 9x = —r ( : ?9t, (49) 

2nnv F2 2-nnvpi -nn(v F i+v F2 ) 

' 1 ( 9l #2 ±9x \ /, n x 

9t= - TTT + i 9t, 50 

2^71 \v F i v F2 v F i+v F2 J 
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where the derivatives are taken with respect to £ = \n(Do/D). By introducing the 
dimensionless coupling constants 



= 1 fm_ + _g2 4ff x . \ y t = ^- ( VF1 + VF2 ^ 2 + 4 ^ F1 ^ 2 (51) 

the four RG equations can be combined into two equations [73] 

v = yt, y't = yyt- (52) 

The flow diagram corresponding to these equations is shown in Figure [TBTb). 

The coupling constant g t describes a combination of two processes: the particles 
transferred between subbands may retain their direction of motion, corresponding to 
momentum transfer ±(k F i — k F2 ), or they may change their direction of motion, 
corresponding to momentum transfer ±(kpi + k F2 ). Both processes are depicted in 
Figure fTBTa). The resulting coupling constant, thus, is proportional to V^{kpi — k F2 ) — 
Vi£(kpi + k F2 ) as given in equation ( 1461) . Consequently, gf' oc k F2 as the density in 
the second subband goes to zero, and one concludes that yf' oc y/v F2 is much smaller 
that y(°\ Therefore the presence or absence of a gap just above the transition to a 
quasi-one-dimensional state is determined by the sign of y^°\ If y^ < 0, the coupling 
constant y t flows to zero, and the system remains gapless. On the other hand, if y^ > 0, 
the coupling constant y t flows to infinity, and the system acquires a gap. 

The interaction constant can be evaluated assuming a Coulomb interaction 
screened by a gate at a distance d much larger than the effective width of the wire, i.e., 

V int (r) = — [-- - 1 | (53) 

1 ' 2e \ t ^r 2 + (2d) 2 J 1 ' 

which is the two-dimensional version of Eq. ([I]) . 

One finds that gi ~ g x ~ 2(e 2 /e) In k F \d (with logarithmic accuracy). On the other 
hand, at densities n 2 <^ 1/d, 

e 2 1 

g 2 (k F2 d) 2 \n- — -, (54) 

e k F2 d 
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Figure 14. Phase diagram of spinless interacting electrons in a quantum wire [65] . 
At any interaction strength there is a finite window of densities where the quasi-one- 
dimensional system supports only one gapless excitation. 

i.e., the interaction constant g^ vanishes in the limit n 2 —>■ 0. This is a manifestation 
of the Pauli exclusion principle: Once the average distance between particles exceeds 
the distance to the gate, the interactions become effectively local. However, identical 
fermions do not interact via a local interaction, hence #2^0. 

Using the above expressions for the interaction constants, one finds the initial value 
y(°> « 3gi/(2nhvFi) > 0. As y^ is positive, the system flows to strong coupling. Thus, 
the system develops a gap close to the transition from a one- dimensional to a quasi-one- 
dimensional state, and therefore the phase diagram Figure [TDT a) does not describe the 
system. 

A second gapless excitation mode appears only once the density is increased further 
beyond the transition point. As the density increases, g2/vp 2 increases and becomes 
comparable to gi/vpx- Then y^ changes sign and eventually one crosses into the regime 
where yt scales to zero. At weak interactions this happens at kp2 ~ 1/(&fi<^ 2 )- Thus, 
at any interaction strength there is a finite window of densities in which the system is 
in the quasi-one-dimensional state but supports only one gapless excitation mode. The 
resulting phase diagram [65] is shown in Figure [TH 

Having found that the behavior at weak and strong interactions is very similar, 
there is no reason to expect that at intermediate interactions no gap exists. One notices, 
however, that the magnitude of the gap strongly depends on interaction strength. In 
the Wigner crystal we find a large Ising gap A ~ \8fi\. At weak interactions, on the 
other hand, the gap scales with a large exponent, A ~ (5fi) a , where a = 1/ (4y(°>). 
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5.3. Intermediate interactions 

The method of choice to treat intermediate interactions in one-dimensional and quasi- 
one-dimensional systems is bosonization. However, bosonization requires a linear 
spectrum. In the present case this is not straightforward because, to describe the 
transition, one is necessarily interested in what happens at the bottom of the second 
subband where a linearization is not justified. 

Alternatively one may bosonize the first subband which has a large Fermi energy 
and keep a fermionic description in the second subband [71] . Thus, the electrons in the 
first subband are described by the bosonic fields <pi(x) and 9i(x) whereas the electrons in 
the second subband are described by the fermionic creation and annihilation operators, 
tp2 ( x ) an d ip2( x )- in particular, 

H = ^I dX ( ( ^ l)2 + ^ (90l)2 ) " £ / dX ^ 2 

- - V(90 1 )(g)V r 12 (g)n 2 (-g) + 7t / alx Ui(>W 2 - c^K^ 1 + h.c.l , (55) 

where ~ e 2 /e. Furthermore, K = (1 + gi/nhvpi) -1 ^ 2 is the Luttinger parameter in 
the first subband. 

The dominant interaction between the bosons and the fermions is the inter- 
subband forward scattering Vu- This coupling can be eliminated by applying a unitary 
transformation 

\K 2 f 1 

■— / dx dy 9\{x)V{x — y)n 2 {y) . (56) 

The new Hamiltonian in terms of the transformed bosonic (<f> c , 8 C ) and fermionic (ipl, ip s ) 
fields then reads 

H v = UHU^ = ^ Jdx (W) 2 + ^m) 2 ^) ~^jdx ^ s 

MM - d^t)e- 2i ^ x) + h.c.l , (57) 



U = exp 



+ lt J dx 

i.e., the inter-subband forward scattering disappears. Comparing the Hamiltonians fl55l) 
and ( 157|) . note that the exponent in the boson-fermion interaction term changes from 
2i9 c to 2\k9 c , where k — 1 — K 2 g x / (irhvpi) — K 2 . It is essential to realize that other 
than that the Hamiltonian preserves its form after the unitary transformation. Namely 
we are still dealing with a plasmon mode coupled to noninteracting fermions. Not only 
the bare interaction in the second subband, but also the effective interaction generated 
by the inter-subband forward scattering vanishes in the limit n 2 — > 0. Additional terms 
that are generated by the unitary transformation can be shown to be irrelevant |74j . 

Since the fermions remain noninteracting, as a next step, they can be bosonized. 
The purely bosonic Hamiltonian could then in principle be subject to an RG approach. 
Or, as it is safe to assume that at intermediate interactions the second mode is still 
gapped, one may use a variational approach instead. One finds that the gap exponent 
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Figure 15. (a) Tunneling into a zigzag crystal, (b) In addition to exciting the plasmon 
mode, the tunneling electron creates a defect. The finite energy cost associated with 
the process manifests itself in a gap in the tunneling density of states. 



decreases with increasing interaction strength until the variational approach is no 
longer valid because the relevant energy scale exceeds the Fermi energy in the second 
subband [66]. For stronger interactions, note that the Hamiltonian (1571) in the limit 
K — > takes the same form as the Hamiltonian of the Wigner crystal with a gapless 
plasmon mode decoupled from the gapped Ising fermionsjjjl 



5.4- Experimentally observable consequences 

As mentioned earlier the computation of observables such as the conductance is 
complicated due to the importance of the coupling to leads. One may speculate, however, 
how the above findings affect observables. 

The experimentally most relevant observable is the conductance. For noninteracting 
electrons, the second subband opens a new channel in the wire and, therefore, at the 
transition the conductance doubles from G = e 2 /h to G = 2e 2 /h (for spinless electrons). 
In the interacting case, however, the second mode is gapped. One might argue [76] that 
the total charge mode (plasmon) remains gapless and, therefore, one should still expect 
a doubling of the conductance at the transition. However, as discussed in section |2l the 
conductance is not determined by the total charge mode only. As the wire is coupled 
to leads, mixing between different channels occurs and, therefore, modes other than the 
total charge mode do affect the conductance. We expect that in this case, too, the fact 
that the second mode is gapped leads to a suppression of the conductance which should 
remain at its one- dimensional value of G = e 2 /h until the second mode becomes gapless 
at a higher density. This means that the transition from a one- dimensional state to 
a quasi-one- dimensional state and the step in conductance no longer coincide. Only 
at higher temperatures T > A does the gapped mode open for transport. Thus, the 
presence of a gapped mode is expected to lead to non-trivial temperature dependence 
of the conductance. 

The gapped mode above the transition should manifest itself most clearly in the 
tunneling density of states. Consider the Wigner crystal limit. In the one- dimensional 
case, the addition of an electron to the system requires the excitation of the plasmon 

% A more careful treatment [75] shows that at strong interactions the weak coupling of the two modes 
is marginally irrelevant and leads to relatively insignificant corrections to the Ising picture of the 
transition. 
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mode in order to adjust the density along the wire. Due to the stiffness of the plasmon 
mode, the tunneling density of states is suppressed: As discussed in section 12.11 the 
Wigner crystal described as an elastic medium can be viewed as a Luttinger liquid, and 
the tunneling density of states of a Luttinger liquid is well known to display a power- 
law suppression at the Fermi level [77J [78]. In the zigzag crystal, the addition of an 
electron to the system also requires to adjust the density along the wire by exciting 
plasmons which suggests a power-law suppression of the density of states. Apart from 
that, however, the addition of an electron creates a defect in the zigzag structure: the 
electron is added to one of the two rows and, thus, interrupts the zigzag pattern as 
depicted in Figure [T5l The energy of such a defect is finite, and, therefore, the density 
of states acquires a gap. 

This behavior is not limited to the Wigner crystal. As discussed in Ref. [76] , at any 
interaction strength tunneling of a single electron into the bulk of the wire excites both 
the gapless and the gapped mode0 and the finite energy cost associated with excitation 
of the gapped mode entails a gap in the tunneling density of states. Consequently 
the observation of a gap opening in the tunneling density of states would allow one to 
identify the transition to the quasi-one- dimensional state. 

6. Conclusion 

The Luttinger liquid physics of one-dimensional electron systems with weak to moderate 
interactions has been studied extensively. The present review focuses on novel 
phenomena due to strong interactions which lead to the formation of a Wigner 
crystal [27]. The strongly interacting regime can be realized experimentally, and 
evidence for Wigner crystal physics has been seen in the conductance of quantum 
wires [791 EQ], the Coulomb blockade peaks in carbon nanotubes [8T] . and possibly [82} 
[83] in the localization features in double quantum wires [84, 85J. While no phase 
transition takes place, at stronger interactions the system properties change due to 
the presence of two very different energy scales, namely the Fermi energy Ep and the 
spin exchange energy J <C Ep. In the strictly one- dimensional regime, one of the main 
features of Luttinger liquid physics is spin-charge separation |7j . In an inhomogeneous 
Wigner- crystal wire, however, spin physics is found to affect the conductance, reducing 
it from 2e 2 /h at T J to e 2 /h in the temperature regime J < T < Ep as discussed 
in Sec. [2j 

Real systems are never strictly one-dimensional, but confined by an external 
potential. The presence of transverse degrees of freedom leads to a transition from a one- 

+ Within the formulation presented here, this can be understood by going back to the unitary 
transformation ([55]) which relates the original degrees of freedom described by Hamiltonian ([55)) 
to the new degrees of freedom described by Hamiltonian (|57[) . In particular, applying the unitary 
transformation to the single electron creation operators ijA, where j = 1,2 is the subband index, one 
may verify that ?/?} as well as contain contributions from both the gapless and the gapped mode of 
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dimensional to a zigzag Wigner crystal at a finite electron density, see Sec. El In contrast 
to the one-dimensional crystal, the zigzag Wigner crystal displays a variety of spin 
ground states as a function of density, see Sec. HI In particular, a ferromagnetic ground 
state — which has been suggested as a possible cause of the conductance anomalies 
observed in quantum wires [T2] — can be realized. While for noninteracting electrons the 
transition to a quasi-one- dimensional state entails the emergence of a second gapless 
excitation mode, this is not the case in the presence of interactions. As discussed 
in Sec. El the orbital degrees of freedom are strongly affected by interactions which, 
for example, are expected to lead to a gap in the tunneling density of states. The 
most interesting properties of the quasi-one-dimensional state in quantum wires are 
summarized in the phase diagrams Figure M (spin properties) and Figure [TH (orbital 
properties). 

Figure M shows the spin phases of the zigzag Wigner crystal obtained under the 
assumption that spin and orbital properties can be treated separately. This approach is 
justified at strong interactions, when the energy scales for spin and charge excitations 
are very different, J <C Ep. However, as interactions become weaker the crystal starts 
to melt, leading to the coupling of spin and orbital degrees of freedom. To study the 
behavior of the system in this regime, one needs to develop a theory that treats spin and 
orbital degrees of freedom on equal footing. This entails a number of open questions: 
Are there remnants of the Wigner crystal phase dominated by the four-particle ring 
exchange in the weakly interacting quasi-one-dimensional state? Is the spectral gap 
discussed in Sec. El robust to the inclusion of spin? 

Figure [TH summarizes the orbital properties of the electron system in a quantum 
wire near the transition from a one-dimensional to a quasi-one-dimensional state. The 
upper line indicates the vanishing of the gap in the second excitation mode. While the 
appearance of a second gapless excitation at a finite distance above the transition has 
been shown in the limit of weak interactions [65], more careful treatment is required 
to explore this phenomenon at finite interaction strength. In the opposite limit of 
strong interactions, with increasing electron density the zigzag regime eventually breaks 
down, giving way to structures with more than two rows. Numerical study of the 
quasi-one-dimensional Wigner crystal [32] shows that the number of rows changes as 
2 — > 4 — > 3 — > 4 — > 5 . . . as a function of density. Since in a wide channel the 
electron density is not uniform across the wire [86l [87] , this trend cannot persist up to 
an arbitrarily high number of rows. Instead one expects that structures with defects 
will have a lower energy. The presence of such potentially mobile defects will be crucial 
for understanding transport properties of quantum wires in that regime. 
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